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Abstract 

Using non-relativistic many body quantum field theory, a master equation is derived for the 
reduced density matrix of a dilute gas of massive particles undergoing scattering interactions with 
an environment of light particles. The dynamical variable that naturally decoheres (the pointer 
basis) is essentially the local number density of the dilute gas. Earlier master equations for this sort 
of system (such as that derived by Joos and Zeh) are recovered on restricting to the one-particle 
sector for the distinguished system. The derivation shows explicitly that the scattering environment 
stores information about the system by "measuring" the number density. This therefore provides 
an important example of the general connection between decoherence and records indicated by 
the decoherent histories approach to quantum theory. It also brings the master equation for this 
system into a form emphasizing the role of local densities, which is relevant to current work on 
deriving hydrodynamic equations from quantum theory. 

PACS numbers: 03.65.-w, 03.65.Yz, 03.65.Ta, 05.70.Ln 
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I. INTRODUCTION 



The notion of decoherence plays a central role in studies of emergent classicality and 
the foundations of quantum theory generally [1]. While it is usually regarded as signifying 
the loss of quantum coherence for a system of interest, it may also usefully be regarded 
as an indication of the degree to which information about the system is stored somewhere 
in the system or in its immediate environment [2, 3]. It is in this way that decoherence 
is related to "generalized measurements". An important application of these ideas is in 
quantum cosmology [4]. There, in applying quantum theory to the very early universe, 
there are no actual measuring device to measure what was happening. The process of 
decoherence, however, guarantees that measurements we make in the present are correlated 
with alternatives in the past. 

These ideas are perhaps most transparent when formulated in terms of the decoherent 
histories approach to quantum theory [2, 4-8]. Other approaches to decoherence, such as 
Zurek's "einselection" approach [1, 9, 10], related density matrix approaches [11] or quantum 
state diffusion [12, 13], may be equally useful for analyzing these issues, but will not be 
explored here. It is the aim of this paper, continuing in the spirit of Ref . [3] , to investigate 
the connection between decoherence and information storage. To fix ideas, we briefly review 
the decoherent histories approach (although the general results of this paper are by no means 
specific that approach). 

In the decoherent histories approach [2, 4-8], probabilities are assigned to histories via 
the formula, 

p(a 1 ,a 2 ,.--)=Tr(C ss pCQ (1) 
where Cg_ denotes a time-ordered string of projectors at times t\ ■ • -t n , 

C a = P^e-^-^P^ ■ ■ ■ e-W^) Pai (2) 

and a denotes the string eti, a 2 , ■ ■ -a n . We are interested in sets of histories which satisfy 
the condition of decoherence, which is that decoherence functional 

D(a,d) = ^{C«pCi) (3) 

is zero when a ^ a'. Decoherence implies the weaker condition that ReD(a, a') = for 
a 7^ a', and this is equivalent to the requirement that the above probabilities satisfy the 
probability sum rules. 
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The stronger condition of decoherence is the more interesting one since it is related to the 
existence of records - information storage about the histories somewhere in the system. More 
precisely, if the initial state is pure, decoherence means that there exist a set of alternatives 
at the final time t n which are perfectly correlated with the alternatives a± ■ ■ ■ a n at times 
ti ■ • • t n [5, 14]. This follows because, with a pure initial state the decoherence condition 
implies that the states Cg\^) are an orthogonal set. It is therefore possible to introduce a 
projection operator Rp (which is generally not unique) such that 



It follows that the extended histories characterized by the chain RpCg^f) are decoherent, 
and one can assign a probability to the histories a and the records (3, given by 



This probability is then zero unless at = flk for all k, in which case it is equal to the original 
probability p(a±, ■ ■ -a n ). Hence either the a's or the /3's can be completely summed out of 
Eq.(5) without changing the probability, so the probability for the histories can be entirely 
replaced by the probability for the records at a fixed moment of time at the end of the 
history: 



Conversely, the existence of records Pi, ■ ■ ■ (3 n at some final time perfectly correlated with 
earlier alternatives a±, • • -a n at ti, ■ • - t n implies decoherence of the histories. 

These issues are most usefully investigated in the context of particular models, and it 
then becomes possible to ask some more precise questions: Which dynamical variables in the 
environment store the information about the decoherent histories? Or what is essentially 
the same thing, how are the "pointer basis" variables stored in the environment? How is 
the amount of decoherence related to the amount of information stored? 

Ref.[3] investigated these questions in the context of the quantum Brownian model 
(QBM), which consists of a particle of large mass M moving in a potential V(x) and linearly 
coupled to a bath of harmonic oscillators. The total system is therefore described by the 
action, 



(4) 






(6) 
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In the traditional discussion of decoherence in this model, it is shown that for a contin- 
uum of oscillators in a thermal state, the influence functional or density matrix become 
approximately diagonalized in position. This may be seen, for example, through the master 
equation for the reduced density matrix p(x,y) of the distinguished system [15], which in 
the high temperature limit is 

In Ref.[3], the issue of how the information about position is stored in the environment was 
addressed. The system is linear in the oscillators, so the classical and quantum dynamics 
coincide for the environment. Classically, the equations of motion of the environment of 
oscillators are, 

m n q n + m n u 2 n q n = -c n x(t) (9) 
The solution to this equation, with fixed p n (0), g n (0) is 

< \ <n\ , • 
Qn{r) = q n (0) cosu n T H smw n T 

m n u n 

f T dtx(t) snWr-i) (10) 

m n u n Jo 

Pn(r) = p n (0) coster - m n u n q n (0) sina; n r 

— C n / dt x(t) COSUJ n (T — t) (11) 

Jo 

where p n = mq n . From this solution, one can see that at the final time r, the positions and 
momenta of the environment of oscillators depend on the particle's trajectory x(t) via the 
temporally non-local quantities 



X* = / dt x(t) smuj n (T-t) (12) 
Jo 

X c n = [ dt x(t) cosLu n (r-t) (13) 
Jo 

These are essentially the Fourier modes of the particle's trajectory x(t). Hence, each oscil- 
lator stores a single Fourier mode of the trajectory, and therefore by using a large number 
of oscillators, information about many Fourier modes is stored from which the approximate 
trajectory may be recovered. Furthermore, since it is the Fourier modes that are naturally 
registered in the environment, rather than positions at each moment of time, decoherence is 



in fact most clearly seen in terms of the variables (12), (13), rather than position, as shown 
in Ref.[3]. The variables are non-local in time so it can only be seen at the level of an influ- 
ence functional expressed in path integral language, rather than a master equation. Hence, 
in this model, it was possible to see exactly how the environment stored the information 
about the system's trajectory in configuration space. Furthermore, a detailed quantitative 
estimate of the information content was also carried out in Ref . [3] . 

Although a very illustrative model, the quantum Brownian motion model is not the most 
relevant model for decoherence in physically interesting situations. Far more physically 
significant is the case in which the environment is a set of light particles which interact 
with the distinguished particle by a scattering process. The resulting master equation, first 
derived by Joos and Zeh, is very similar in form to the QBM case, Eq.(8) [11, 16, 17]. But 
the dynamics of the environment, and therefore the means of information storage, are rather 
different. 

The aim of this paper is to investigate the connection between decoherence and records 
in the case of decoherence by a scattering environment. In some ways it is simpler, since, in 
the usual assumption of widely separated timescales for system and environment dynamics, 
each environmental particle scatters briefly off the distinguished particles, and moves freely 
thereafter, carrying some information about the distinguished particles. This process can 
therefore be described by a Markovian master equation, and the process of information 
storage and decoherence may be described in a moment by moment manner (unlike the 
quantum Brownian motion case, where the environment oscillators store information about 
the entire history). This means in fact that we do not need to make use of the full machinery 
of the decoherence functional - it is sufficient in fact to look at the evolution of the reduced 
density operator. 

In the quantum Brownian motion model case, the system variables the environment 
measures were actually identified quite simply, from the classical equations of motion. In 
the scattering case, the system variables measured most directly by the environment are also 
determined quite easily, by examining simple scattering processes. In particular, suppose 
we consider the scattering of some light particles off a dilute gas of a set of more massive 
particles with coordinates q^. Then it follows quite straightforwardly from simple scattering 
theory (and we will in fact demonstrate this) that the scattering amplitude is proportional 
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to the Fourier transform of the number density of the massive particles, 



(14) 
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This means that, loosely speaking, for a known interaction potential, measurements of the 
initial and final momenta of the scattering environment determine the number density of 
the distinguished system. 

Of course, the number density is closely related to position, which is normally held to 
be the preferred basis in these calculations. But, following the lead of the oscillator model, 
we expect decoherence to look simplest in terms of the dynamical variables which are most 
simply and directly stored in the environment. Our aim is therefore to give a derivation of 
the master equation which emphasize the central role played by the number density. We 
have found that the derivation is in fact most transparent in terms of non-relativistic many 
body quantum field theory, where the number density appears very naturally. We will give 
an alternative and more general derivation of the master equation, using many body theory, 
which brings out the role of local number density more clearly, hence showing the connection 
with records. 

It is pertinent at this stage to mention the Lindblad form of the master equation [18], 
which is the most general possible form a master equation can take under the assumption 
that the evolution is Markovian (a condition well-satisfied in a wide variety of interesting 
models). The Lindblad master equation is 



Here, H is the Hamiltonian of the distinguished subsystem (sometimes modified by terms 
depending on the Lj) and the n operators Lj model the effects of the environment. The 
master equation of quantum Brownian motion, for example, is of this form with 



as described in Refs. [13, 19]. (Actually, the master equation (8) is not strictly of the Lindblad 
form, and as a consequence can suffer from a violation of positivity [20]. However, the 
difference between Eq.(8) and the Lindblad form with L given by Eq.(16) is of order 1/T 
which does not matter for high temperatures). 




(15) 




(16) 
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The Lindblad operators Lj determine the basis in which the density operator tends to be- 
come approximately diagonal, or what is essentially the same, the sets of variables describing 
an approximately decoherent set of histories. This may be seen from the formal solution to 
the Lindblad equation [21]. Consider the case of a single Lindblad operator L = Lr + iLj, 
where Lr, Lj are hermitian. Divide the finite time interval [0, t] into K subintervals, so 
that t = KSt, and let St — > 0, K — > oo, holding t constant. Then, the formal solution to 
the Lindblad equation is obtained by taking the limit St — > 0, K — > oo (with t fixed) of the 
expression, 

K 



Pit) = Jd%...d 2 l K 



x fl exp t^(t m L - £ m L^)\ expf-y \L-£ m A exp (~H'6t ) p(0) 



m=l 

x n ex p (i H ' 5t ) ex p (-| i l - ^i 2 ) ex p (- y (c^ - w f )j (17) 

Here, H' = H + ^[L, L^j, and the £ m are complex numbers at the discrete moments of time 
labelled by m. We use the notation 

\L -C| 2 = - ReC) 2 + (i/ - ImC) 2 (18) 

The ordering of the operators at each moment of time is irrelevant in the limit St — > 
(although the operators at different times are time-ordered, according to increasing to). 
That this is the solution is readily verified by explicit computation. The solution has the 
form of a "measurement process" of the L's, continuous in time, with "feedback" via the 
terms {t m L-£ m V) [22]. That is, one can see that the effect of the environment is to produce 
a tendency towards diagonality in L. 

We shall show that a many-body theory derivation of the master equation in the case 
of a scattering environment leads to a master equation of the Lindblad form (under the 
assumption that the environment dynamics are much faster than the system dynamics) , and 
that the Lindblad operators are essentially the local number density. The previous forms 
of the master equation are recovered in the one-particle sector for the system of massive 
particles. 

This work grew out of a more ambitious programme, in the context of the decoherent 
histories approach to quantum theory, which aims to give a very general account of emergent 
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classicality. In particular, it is asserted that at sufficiently coarse-grained scales, the local 
densities (number, momentum, energy) define a a set of habitually decohering variables, 
even in the absence of an environment [2, 23]. This is because they are locally conserved, 
and therefore slowly varying when coarse-grained over sufficiently large volumes, and thus 
are expected to be approximately decoherent (because exactly conserved quantities are ex- 
actly decoherent in the histories approach). This is therefore a different mechanism for 
decoherence than the usual mechanism of decoherence through an environment. Hence, in 
order to close the gap between the familiar system-environment models and the less fa- 
miliar hydrodynamic models without an obvious environment, it is useful to rewrite the 
system-environment models in terms of local densities as we do here. 

In Section 2, we briefly review many body field theory, and carry out the simple scattering 
calculation leading to the result the the scattering particles effectively measure the local 
number density, Eq.(14). 

In Section 3, we use many body field theory to derive the master equation for the system, 
using a slow motion limit for the gas of massive particles. As anticipated it has the Lindblad 
form in with the Lindblad operators proportional to the Fourier-transformed number density 
N k . 

In Section 4 we show that our master equation reduces to an earlier result of Gallis and 
Fleming [17] in the one-particle sector for the gas of massive particles. (This is essentially 
the same as the master equation of Joos and Zeh [11] but the comparison with Gallis and 
Fleming is more direct). 

The master equation of Sections 3 and 4 does not have any dissipation and is analogous 
to the Lindblad equation of quantum Brownian motion with L proportional to x. In Section 
5, we go beyond the slow-motion limit to derive a master equation with dissipative terms. 

We summarize and conclude in Section 6. 

II. MANY BODY FIELD THEORY 

The dynamics of a many body system is very conveniently handled using many body 
quantum field theory. We now set up the formalism of many body field theory [24, 25] 
which we will use to derive the master equation. We consider a set of non-relativistic system 
particles described by a field if)(x.) interacting through a potential </>(x) with an environment 
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described by a field x(x). The total system is described by the Hamiltonian 
H = jd 3 x (^V^(x) • VV(x) + lvx f (x) • V*(x)) 

+ l -J d 3 xd 3 x' ^t( x )^( x ')0( x _ X ')x t (x , )x(x) (19) 

(For simplicity we set ft = 1 hereafter). In this language, the number densities iV(x) and 
n(x) of the system and environment fields are 

iV(x) = ^ f (x)^(x) (20) 
n(x) = x f (x)x(x) (21) 

(This is the field-theoretic version of Eq.(14)). 

The above relations are also more conveniently written in terms of a k and fe k , the anni- 
hilation operators for the system and environment, respectively, and the Hamiltonian then 
is 

H = E { E A a « + ( 22 ) 
q 

+ ^77 E ^( k 2 -k 2 )4 1 C a k' 1 ^ 

k'j+k^ki+ka 

where _E q = q 2 /2M, u; q = q 2 /2m, V is the spatial volume of the system (which we assume 
is in a box) and 

i/(k) = y d 3 x e" ik - x 0(x) (23) 
The Fourier transformed number densities are 

iV k = E a W+k (24) 
q 

«k = E & q & q+k (25) 

q 

and one may see that the Hamiltonian has the more concise form 



H = E ( E A a « + "q&W + ^7 E *(k)iV k n_ k (26) 
= H + H mt (27) 



From these relations we see that the environment couples to the number density of the 
system. It is this feature of many body field theory that makes it the appropriate medium 
for the derivation of the master equation emphasizing the role of number density. 
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The S-matrix is 



S = T cxp (-i |°° dt H mt (t^j (28) 



where 

^(^k^W^W) ( 29 ) 



2^ k 



and here 



Mt) = E a q« q+ ke 1 ^-^)' (30) 
q 

«k(*) = E fe iwe !K ^ +k)t (31) 

q 

We may now use this formalism to look at a simple scattering situation to determine how 
the environment stores information about the system. In the quantum Brownian motion 
case, the nature of information storage was determined in essence by solving the classical 
equations of motion. A similar strategy works here too. Let us suppose the distinguished 
system is classical, and consider what happens when the environment scatters off it. Suppose 
the environment starts in an initial momentum state |ko) and scatters into a final state |k/). 
The scattering amplitude for this process, to first order, is 

j roc 

(k/|S|k > = — j_Jt^(k)N k (t){k f \n_ k (t)\ko) 

/(k) J dt N k (t) e ^ k /-^o)* (32) 
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where k = kj — k . This simple result shows that a single scattering event by the environment 
stores information about the Fourier transform (in space and time) of the number density. 
It is in this sense that the number density has a preferred status - this is the variable that 
is measured most directly by the environment and is the exact analogy of the relations 
Eqs.(12), (13) in the quantum Brownian case. 

The measured variables above are of course non-local in time, involving a temporal Fourier 
transform of the number density, so cannot in fact be compatible with a Markovian master 
equation. Under a reasonable slow motion assumption, the system timescale is much slower 
than the environment timescale, and we may ignore the time-dependence in N k (t), yielding 

(k/ISIkb) = ^ i/(k) N k 5(cu kj - cu k{) ) (33) 

This corresponds more directly to a Markovian master equation, as we shall see. 
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It remains to briefly sketch the connection between these results and the discussion in 
Section 1 of records in the decoherent histories approach. We imagine that the environment 
consists of a very large number of particles which scatter off the system particles. Each 
scattering event consists of an incoming environment particle with momentum k scattered 
into a final state of momentum kf, as outlined above. After the scattering event, which 
is essentially instantaneous (on the timescale of system dynamics), we may imagine that 
the scattered particle propagates freely and may be measured at any time in the future. 
Therefore, the records in the decoherent histories approach consist of projections at the end 
of this histories onto the momenta of all the scattered environment particles, from which 
the number densities of the system at a series of times may be retrodicted. 

III. DERIVATION OF THE MASTER EQUATION 

Following the method first used by Joos and Zeh [11], we may derive the master equation 
for the reduced density operator p of the system by considering the scattering of the envi- 
ronment off the system, to second order in interactions. We assume that the system and 
environment are initally uncorrelated, so the total density operator is 



We also assume that each scattering event takes place on a timescale which is extremely 
short compared to the timescale of system dynamics. This means that in an interval of time 
At which is long for the environment but short for the system, we may write, 



(where we are using the interaction picture). Expanding (35) to second order, the S'-matrix 
may be written, 



Pt = Po ® Pe 



(34) 



p T (t + At)=Sp T (t)S^ 



(35) 



S = l + iU 1 -U 2 



(36) 



where 



/. 



oc 



u 1 = - 



dt H int (t) 



(37) 



— oo 



and 




(38) 
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The requirement of unitarity, S 1 = implies that Ui = U\ and 

U 2 + U\ = U\ (39) 



We will therefore write 



where B = B\ so we now have 



U 2 = hjl + iB (40) 



S=l + t{U 1 -B)-hjl (41) 

Inserting this in (35), we obtain 

^-At = i[U 1 - B, p T ] + \J lPT U x - l -U 2 lPT - l -p T U 2 x (42) 

We now trace Eq.(42) over the environment to obtain the master equation for the system 
density operator p. As is usual in this sort of model, we assume that the environment is so 
large that its state is essentially unaffected by the interaction with the system. Since the 
total density operator starts out in the factored state (34), this then means that, to a good 
approximation, px persists in the approximately factored form p £g> ps, and we may insert 
this in the right-hand side of Eq.(42) [26]. We thus obtain the preliminary form for the 
master equation 

^ At = i[rr e (U lP e) - Tr g (Bp e ),p] + Tr £ (u^U^ - \u 2 iPt - \p T u() (43) 

We now work out these terms in more detail. We first consider the simple but useful slow 
motion approximation, in which we ignore the time-dependence of N k {t). This implies that 

~ £ Kk)iV k J2 & q & q-k 2vr5K - cu q _ k ) (44) 

ZV k q 

The important terms for decoherence are the final three terms on the right-hand side of (44). 
When traced, these give, 

Tr £ (ViPTfA - l -U 2 lPT - \ptUI) = £ c ( k > k ') (^k'piVk - \n v N vP - ^pN k N k ,^j (45) 

where 

C(k, k') = U(k)v(k!) £ <K^q - ^q-k)5(cJ q / - Wq'-k') (& q &q-k& q ' &q'-k')f ( 4 6) 

qq' 
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We will take the environment to be a thermal state, which is diagonal in the momentum 
states. It follows that 

(&] 1 &q-k&q'&q'-k')£ OC 5 q , q '_ k ' ^q',q-k (47) 

This implies k = — k', and also that the two delta-functions are the same in Eq.(46). We 
then interpret the square of the delta-function in the usual way, 



[5(cjq - CJ q _ k )] 2 = 5(0) <5(Wq-^ q _ k ) 
At 

~ 2^ 



At 

— 5(Wq-Wq_ k ) (48) 



We now have 



where 



c(k,k') = 4,- k 'c(k)^ (49) 



c ( k ) = ^K k )l 2 ]C ^K-^kK^Vktk^f 

= ^IKk)| 2 E 5K-c, q _ k )(6 q 6 q ),((6t_ k&q _ k ), + i) (50) 

The terms involving environment averages have the usual thermal form (for a bosonic envi- 
ronment), 

<HN>« - ^ _ i < 81 ) 

where f3 = 1/kT with T temperature, and /x is the chemical potential. 

The form Eq.(49) means that the important terms in the master equation are of the 
Lindblad form, 

Tr^ (u^tU, - 1 -U 2 iPt - \ptV\) = At E c ( k ) ( N *P N t - ~ \pN^ (52) 

where we have used the fact that = AL k . The remaining two terms in Eq.(43) clearly- 
just modify the unitary dynamics of the system. First we have 

Tr^ (U lP£ ) = i- E Kk)iV k E(^-k)^ 2tt5K - cu q _ k ) (53) 

AV k q 

Clearly from the term (& q & q _ k )£ this expression will be zero unless k = 0, and therefore it is 
proportional to N, the total particle number operator (although the overall coefficient will 
need to be regularized). This therefore contributes a term to the master equation of the 
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form [N, p] . We assume that there is a fixed number of system particles so it is reasonable 
to take this term to be zero. 

The other remaining term in Eq.(43) involves the time ordering terms in £7 2 and is a 
bit more complicated to evaluate. Fortunately, the detailed form of this expression is not 
needed here, and it can in fact be easily shown that this term has the form 

Tr e (B Pe ) = At£ d(k) A k iV k (54) 

k 

for some coefficient d(k) which we will not need. Inserting all these results in Eq.(43), the 
factors of At all drop out, and we obtain, in the Schrodinger picture, 

^ = -*[#o-E d(k) iV k ATt, p ] + £ c (k) [N kP Nl - ^NlN kP - \pNlNC) (55) 

As desired, this is the Lindblad form with the Lindblad operators given by 

L k = c5(k)iV k (56) 

We have therefore produced a derivation of the master equation for a scattering environ- 
ment which shows very clearly the connection between the preferred basis (diagonalization 
in the Lindblad operators) and the information storage about the system, as indicated by 
the simple scattering calculation, Eq.(33). 

It is interesting to note that the decoherence effect is second order in interactions, but 
we were able to anticipate it from the simple first order calculation, Eq.(33). The reason for 
this is the relationship Eq.(40), which shows that the important part of the second order 
terms is the square of the first order terms, and this is a consequence of unitarity. 



IV. COMPARISON WITH PREVIOUS WORKS 



It is useful to check that the master equation we have derived reproduces known results 
when we restrict to the one-particle sector for the system. We will compare the results of 
this to the derivation of Gallis and Fleming [17] (which is essentially the same as Joos and 
Zeh [11] and Diosi [16]). 

In the one-particle sector we may work with a density matrix p(k, k') = (k|p|k'), or 
equivalently p(x, y) in the position representation. We use the relations 

(57) 

[A q ,a k ] = -a k+q (58) 
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A q ,a k 



These relations imply that 

iV q p(k, k')AL q = p(k - q, k' - q) (59) 

iV_ q iV q p(k,k') = P(k,k') (60) 

P (k,k')iV_ q iV q = p(k,k') (61) 

In the position representation, this means 

iV k p(x, y)AU = e tk -( x - y V(x, y) (62) 

The master equation for the one-particle density operator p(x, y) is then 

dp(x, y) 
dt 

where 



= -i(x| [if„, P] |y> - F(x - y)p(x, y) (63) 



F ( x - y) = 2F2 / rf3 ^ 3A; l z/ ( k )! 2 ^K-k + !) 5 K - "q-k) (! " e ik - (x - y) ) (64) 

Note that the term involving the coefficient d(k) in Eq.(55) drops out because [N^N^, p] = 
in the one-particle sector. 

To compare this with the Gallis-Fleming result [17], we first introduce the quantity 

/(k,k') = |^(k-k') (65) 

(which appears in the usual Born approximation to first order scattering). Then, letting 
k — > — k + q in (64), we get 

A 2 

F(r) = -^ / d 3 qd 3 k |/(q,k)| 2 n q (n k + 1) % q - u; k ) (1 - e^'^) (66) 
m z J 

The delta-function implies that q 2 = k 2 , and we find that 

F(r) = ^ f dqq 3 n q {n q + l) J dQdQ'\f(q, k)| 2 (l - e^"*)-) (67) 

This in fact agrees with Gallis and Fleming if we identify q/m as their v(q), the speed of 
the incoming particles, and 4n 3 q 2 n q (n q + 1) as the density of particles with speed q. In the 
one particle sector there is therefore agreement with earlier work. (At least up to an overall 
numerical factor which we could not rectify. However, we have also spotted some small and 
probably insignificant numerical errors in Ref. [17]). 
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Mention should also be made of the master equations derived by Unruh and Zurek, which 
used a field as an environment for a particle [27] , and Anastopoulos and Zoupas [28] , which 
used a photon field and as an environment for a spinor field. Also of relevance is the general 
account of the derivation of master equations given by Omnes [29]. These works are rather 
different to the present paper. 



V. BEYOND THE SLOW MOTION APPROXIMATION 

The derivation above assumed, in essence, that the system dynamics are infinitely slow. 
Not surprisingly, the resulting master equation does not involve dissipation, since, in the 
approximation used, the system is essentially at rest for the timescale of a single scattering 
event. It is analogous to the master equation of quantum Brownian motion with the Lindblad 
operator L proportional to x, Eq.(16). To get a more realistic equation with dissipation we 
therefore need to go beyond the slow motion approximation. 

Because the local number density is a locally conserved quantity, it obeys a continuity 
equation of the form, 

A k = -ik • P k (68) 
where P k is the local momentum density, 

p k = E(q+^ k ) a q a q+k (69) 

It is reasonable to expect that the master equation will involve this operator when we go 
beyond the infinitely slow limit. We now briefly repeat the derivation of the master equation, 
this time allowing a slow time- dependence in N k (t). 
We have 

Ui = ^rEKk) E*&-k / dt N k (t) e^-^X (70) 

ZV k q J 

To take into account the time-dependence of A k (t), we write, 

N k (t) = N k + tN k + --- (71) 

where A k is given in terms of the momentum density, Eq.(69). Inserting this in Eq.(70), the 
factor of t may be rewritten in terms of a delta-function derivative, yielding, 

U i = ^77 E Kk) E & q (A k - cu q _ k ) - *A k 5'K " wq-k) + • • •) (72) 

ZV k q 
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We now use this expression for U\ in the derivation of the master equation. So for example, 
we get, in place of Eq.(45), 

Tr^ (Ulpe) =AtJ2 c(k) N^p + % £ 6(k) (iV k iV k - iV* AT k ) p + • • • (73) 

k k 

where c(k) is given by Eq.(50) and 

6(k) = |z/(k)| 2 £ % q - c q _ k ) 5'K - Wq _ k ) (6 q 6 q ) £ ((6 q _ k 6 q -k) £ + l) (74) 
q 

This coefficient may in fact be shown to be simply related to c(k). The delta- function 
derivative may be dealt with by noting the formal relation 

S(x)^x) = ~[6(x)] 2 (75) 

Now note that 

- cu q _ k = ^ (2k • q - k 2 ) (76) 

It follows that the delta function derivatives may be expressed in terms of derivatives with 
respect to as 

5'(cj q - CJ q _ k ) = -j^ fcj — 5(cj q - CJ q _ k ) (77) 

Inserting these relations in Eq.(74) and integrating by parts yields, 

W = -^(k)| 2 E [% q -c, q _ k )] 2 |^|- (» q >,((t k V^ + l)) (78) 

where we will interpret the square of the delta-function as in Eq.(48). 
Now for simplicity work in the high temperature limit, so 

(&UVk>£ « 1 (79) 



and 



It follows that 



(b^s « e^e~^ (80) 



= -P— (b%} £ (81) 

7TZ 
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and we arrive at the very simple result, 




(82) 



It is not difficult to see that we then arrive at a master equation which is once again of the 
Lindblad form, but this time with Lindblad operators of the form 



(up to terms of order ft 2 , which can be dropped in the approximation we are using). This is 
clearly closely analogous to the QBM result, Eq.(16). (A closely analogous formula appears 
in Diosi's paper [16]). 

VI. SUMMARY AND DISCUSSION 

We have given a derivation of the master equation describing a many-body system inter- 
acting with a reasonably general class of environments. The form of the master equation 
emphasizes the central role of the local number density, which is the system variable mea- 
sured most directly by the environment in a scattering situation. 

We did not in fact give a specific form for the interaction between the system and en- 
vironment since it was not necessary to illustrate the general points we are making. Some 
specific forms for this interaction are discussed elsewhere [11, 17]. 

The derivation reduces to familiar results of Gallis and Fleming [17], Diosi [16], and Joos 
and Zeh [11], when we restrict to the one-particle sector of the many-body field theory. 
The many-body derivation confers some advantages of the usual derivations (which consider 
scattering theory in quantum mechanics) in that they avoid essentially classical assump- 
tions about fluxes of scattering particles. Our derivation also has the possibility of being 
extended to a low temperature regime (and to Bose-Einstein condensation, for example) and 
to fermionic environments, although we do not discuss this here. 




(83) 
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